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Abstract 

I theoretically found eigenmodes and growth rates of relativistic current filamentation instability 
in collisional regimes, deriving a generalized dispersion relation from self-consistent beam-Maxwell 
equations. For symmetrically counterstreaming, fully relativistic electron currents, the collisional 
coupling between electrons and ions creates the unstable modes of growing oscillation and wave, 
which stand out for long-wavelength perturbations. In the stronger collisional regime, the growing 
oscillatory mode tends to be dominant for all wavelengths. In the collisionless limit, those modes 
vanish, while maintaining another purely growing mode that exactly coincides with a standard 
relativistic Weibel mode. It is also shown that the effects of electron-electron collisions and thermal 
spread lower the growth rate of the relativistic Weibel instability. The present mechanisms of 
filamentation dynamics are essential for transport of homogeneous electron beam produced by the 
interaction of high power laser pulses with plasma. 

PACS numbers: 52.25.Fi, 52.27.Ny, 52.35.Qz 
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I. INTRODUCTION 

Relativistic laser-plasma interactions have been a topical issue for the past decade , in 
the context of ignitor physics of inertial confinement fusion, aiming at additional fast heating 
and subsequent ignition of highly compressed targets by means of an external intense laser 
pulse jjj. The laser pulse drives relativistic currents, compensating return currents, and 
creates a pattern of counterpropagating currents which are subject to current filamentation 
instabilities (CFI) including the Weibel mode |3|. In the theoretical arena, numerous ver- 
sions of analytical and numerical methods have been developed in the past, to explore this 



type of instabilities 
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In general, the physical 



mechanism of the electromagnetic CFI is explained as follows: When the compensation of 
the counterpropagating electron currents is disturbed in the transverse direction, magnetic 
repulsion between the two currents reinforces the initial disturbance. As a consequence, a 
larger and larger magnetic field is produced as time increases, degrading the transport prop- 



erties. Many efforts have been devoted to this crucia 



electron beams 
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problem, both related to laboratory 



as well as to astrophysics [2JJ, 2^, 23]. Concerning laser interac- 
tion with plasmas, for the low (nonrelativistic) intensity regimes the inverse bremsstrahlung 
absorption is known to predominate around the cutoff region, where the Weibel-type in- 
stability associated with temperature anisotropy can take place 24I 
plasmas, the collisional and related nonlocal effects were investigated 

The original motivation for this work was triggered by more recent publications that have 



25 1 . In the ablative 



been quantitatively treated with the counterstreaming relativistic CFI in the collisionless 



limit |27, 



28, 



29|. A series of works could be linked with the ignitor physics by irradiating a 



relativistic laser pulse: Fast ignition of the compressed fuel requires at least 10 — 100 kJ of 
external energy to be deposited within ~ 10 ps into the precompressed core 3(|. If carried 
by 1 — 10 MeV electrons, it implies a current of 0.1 — 1 GA which exceeds the transport limit 
of about 100 kA [3jJ, by more than a factor ~ 10 3 . The essential feature is the breakup of the 
relativistic electron beam into many filaments when propagating in dense plasma [3, l^ . 
The physics underlying this phenomenon is just the CFI as mentioned above. This type of 
instability leads to nonlinear filamentation and coalescence of the relativistic electron beam 
flflfl HHQ and the form a ti „ n of strong magnetic ne. ds UU HQ" 



In the 



astronomical point of view, one often encounters such morphology in a variety of celestial 
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objects, particularly, in the astrophysical jets [41]. State-of-the-art observations by utilizin 



very long baseline interferometry have revealed the filamentary structure of the jets 42], 
involving transverse magnetic fields More recently, large-scale toroidal magnetic fields 
have been discovered in the Galactic Center [44j], accompanied with splendid filamentary 



radio arcs 
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Regarding the ignitor physics of laboratory plasmas, the effects of collisions and beam 
thermal spread play a significant role in the CFI caused by the ultrahigh relativistic currents 
through ablative coronal plasma [37], whose density rises from n c m 10 21 cm~ 3 (cutoff 
density) to ~ 10 26 cm -3 (a thousand times solid density), over a radial distance of ~ 10 2 /zm 
The role of collisions in a laser-produced plasma has been discussed in the early 



literature by Motz 48], though the description was restricted to the nonrelativistic fashion. 
It is important to note that the collision frequency invokes additional parameter disturbing 
the universal density scaling, in terms of the time scale of a plasma oscillation period uj~^ 
and the spatioscale of a skin depth c/uo pe , which is valid only for collisionless regimes. In 
fully relativistic regimes, the collision cross section should be evaluated by using the Mott 
scattering formula |4j| . Presuming a small angle scattering and averaging over the angle yield 
the electron-ion collision frequency that can be defined as v ei = (riiY ei /c 3 )(T/ /z 3 ), where rij is 
the number density of ions, Y ei = 47r(Ze 2 /m ) 2 lnA, T = (l + /i 2 ) 1//2 , ii = p e /m c, and p e , m , 
Z, and InA are the electron momentum, the electron rest mass, the averaged charge number, 
and the Coulomb logarithm, respectively 50], and the electron-electron collision frequency 
is given as v ee = (2/Z)u ei . Introducing the current neutral condition of rib/n p ~ \v p /c\ <C 1, 
where rib, n p, and v p are the beam electron density, the plasma electron density, and its 
velocity, respectively, we find the ratio of v ei /uj pe ~ 27r 1 / 2 Z 9 / 2 e 3 ?2^ 2 (m c 2 )~ 3 / 2 (n i /nft) 3 lnA 
for the plasma return current. When assuming rib ~ n c w 10 21 cm -3 to be almost constant, 
the ratio can be estimated as 

u pe V10 23 cm- 3 ; \ 5 ) ' K J 

Equation indicates that indeed, collisional effects are important in the supersolid density 

regions of rii > 10 23 cm -3 , and particularly, for the Lorentz plasma with Z>1. In addition, 
the beam thermal spread will occur there, since the electrons penetrating through the ab- 
lative corona are expected to be thermalized via the collisional and collisionless dissipative 
processes j3]. The thermal effects involved with the ratio of the transverse temperature to 



the total energy of beam electron, i.e., ~ Tf,±/ (Tm c 2 ) with expected values in the range 
less than unity 3?j, which also disturbs the aforementioned scaling. 

In this paper, I present fundamental eigenmode properties and growth rates in the linear 
stage of the CFI including the effects that violate the universal density scaling. The theory 
is expanded in fully relativistic regime on the basis of the self-consistent beam-Maxwell 
equations. To the best of my knowledge, a comprehensive treatment of electron-electron 
and electron-ion collisional effects on relativistic counterstreaming electron currents with 
thermal spread has not been carried out, so far. For the symmetrically counterstreaming 
currents, I found that collisional and thermal effects are likely to lower the growth rate of 
the relativistic Weibel instability. The most significant result is that the finite collisional 
coupling between electron and ion creates the growing oscillatory and the growing wave 
modes, which stand out for long-wavelength perturbations, and in the moderate to strong 
collisional regime, even for short- wavelength perturbations, the growth rate of the oscillatory 
mode exceeds that of the suppressed Weibel instability. In this aspect, the present work goes 
beyond the framework of the well-established theory of electromagnetic instabilities. The 
asymmetric configuration effects of the counterstreaming currents are also investigated by 
using a slow return current approximation. Furthermore, I argue that the collisional coupling 
between electron and electron creates a growing wave mode, but its growth rate is lower than 
that of the suppressed Weibel instability. It is also shown that thermal effects participate 
in lowering the growth rate of the Weibel instability. Although, in the case that includes 
thermal corrections, the present calculation is valid for the smaller wave number k as shown 
later, the most interesting range k ~ u pe / c can be fairly covered. 

In order to spell out these subjects, the present paper is organized as follows: In Sec. HTJ 
the linear theoretical analysis of the relativistic CFI is expanded systematically. The basic 
equations introduced in Sec. Ill Al are linearized so as to obtain a dispersion relation of the 
CFI along the manner outlined in the Appendix. The generic dispersion relation is pre- 
sented in Sec. IIIB| and its approximate expressions are derived in Sec. Ill CI In Sec. for 
an application, the newly derived equations are solved for typical parameters of counter- 
streaming relativistic currents, and the properties of complex eigenmodes are investigated 
for the cases including the effects of electron-electron collision fSec. lIII A"| . electron-ion colli- 
sion fSec. IIII B|) . thermal corrections (Sec. IIIIC^ . and for the special case without collisional 
and thermal effects (Sec. IIIIDj) . as well as, for the case with an asymmetrical configuration 
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of counterstreaming currents (Sec. IIIIE|) . closely relevant to ignitor physics. As a matter of 
convenience, formulas of the growth rates are explicitly written down for some interesting 
cases. Finally, Sec.|TV]is devoted to concluding remarks. 



II. GENERALIZED DISPERSION RELATION OF THE RELATIVISTIC 
CURRENT FILAMENTATION INSTABILITY 

A. Basic equations and assumptions 

Begin with the nonlinear beam-Maxwell equations that include both friction and pres- 
sure terms. Assuming the ions to be at rest and to provide a uniform charge-neutralizing 
background, we study the relativistic dynamics of two uniform, counterstreaming electron 
currents by employing the following set of equations in the dimensionless form, 

^-V.j. = 0, (2) 

dp a V P a 

~^7~ + (v a • Vp a ) = — (E + v a x B) — u ei p a - v ee (p - p a ) -, (3) 

VxE = --, (4) 

VxB = ^ + Ej- ( 5 ) 



V-E=l-^n a , (6) 

a 



where v a = p a / a/ 1 + pi, j a = —n a v a , and the subscript a = 1, 2 labels the two electron com- 
ponents and a labels the countercomponent of a, viz., a — 2, 1 for a = 1, 2, respectively. In 
particular, z/ e j and v ee describe the electron-ion and electron-electron momentum exchanges, 
and the other notations are standard. For normalization, I have used the initially uniform 
density n , the speed of light c, and the electron plasma frequency u pe = ^Annoe 2 /m . 
Note that the Poisson Eq. © is equivalent to a combination of the continuity Eq. (j2J) and 
the Ampere-Maxwell Eq. . 
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According to a procedure similar to that developed by Califano et al. [2jj, I investigate 
the behavior of small amplitude perturbations by linearizing Eqs. — ©• I impose current 
neutrality of Yl a n o,a v o,a — 0, where v 0ja = ^o,aX are the initial velocities in x direction. 
Under the current neutrality, there exists no magnetic field initially. The CFI is studied 
in the x-y plane. In order to derive the dispersion relation, all perturbed quantities are 
assumed to be in the form of T(y,t) = J r ex\i[i{kyy — ut)}. As a result, a magnetic field 
B = (0,0, B z ) and a corresponding electric field E = (E x ,E y ,0) are generated. The purely 
transverse mode is investigated in detail throughout this paper, i.e., k x = and k y = k ^ 0, 
whereas the longitudinal mode such as two-stream instability with k x ^ and k y = is not 
taken into account at the moment. The pressure of each electron components is connected 
with its density by a polytropic relation, which depends on the characteristic frequency uj 
and wave number k of the mode being considered. As it is well known, in the case that 
the ratio of u/k is much larger than the electron thermal speed, the adiabatic exponent of 
7 e = 3 is adequate for the polytrope Henceforth, I assumed P a /n\ = const, so that 

VP a = 3T 0>a Vri a . Then, we get a closed form of the linearized Eqs. (J2J) — ©• The background 
ions are supposed to be fixed on the short time scale of ~ u~ e . Along these assumptions, 
in the nonrelativistic limit Eqs. © — © involve the resistive transverse wave modes and 
longitudinal wave mode with thermal correction which were discussed in Ref. (48 1. 

B. The dispersion relation including collisional effects and thermal corrections 

The extended dispersion relation of the relativistic CFI including collisional and thermal 
effects is then found self-consistently. In collisionless cases the dispersion relation that can 
be expressed as functions of uj 2 and (k/co) 2 does not include the imaginary unit i 



explicitly and contains the purely real and purely imaginary solutions of uj{k) 27[. The 
purely real solutions consist of the pairs of the positive and negative solutions, correspond- 
ing to purely oscillatory and/or purely oscillatory wave modes, while the purely imaginary 
solutions are concomitant with the complex conjugate solutions, to yield purely growing and 
purely decaying (damped) modes. 

In the collisional case considered here, the dispersion relation includes the imaginary unit 
explicitly. Hence, the solutions of uj{k) may depart from the real and imaginary axis in the 
complex plane. In this sense, hereafter we refer to such solutions, i.e., complex eigenmodes 
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with real and imaginary part, as dephasing modes. Below, I explicitly write down the 
generalized dispersion relation containing the dephasing modes. After some manipulations 
outlined in the Appendix, the dispersion relation can be obtained in the complex form of 
Re[u>, (k/uj) 2 } + ilm[uj, (k/ou) 2 } = 0, where 

Re[u, (k/u) 2 } = [(l + rQ k 2 ) - (Q u 2 + Q u 2 T Q k 2 )} 

x [(i + m k 2 ) (i - n k 2 ) - (n 2 2 + n 21 2 T n k 2 ) - (n£ + n 31 2 T n k 2 ) n k 2 } 

— (yi 12 + ^l2,T^fc ) [(^22 + ^22,T^fc ) + (^32 + ^32,T^fc ) ^A; ] 



[(f2 41 + fi 41T fi fc ) (fi 43 + ri A3T Q k ) (fi 42 + {l A2 T {l k ) (f2 44 + ^44^0^ )] £l k , (7a) 

im[^, (k/u) 2 ] = - {q u 2 + n~l T n k 2 ) 
x [(i + rnf) (i - nf) - {n 2 2 + n 21 2 T n k 2 ) - (n 3 2 + n^nf) n k 2 ] 
- [(1 + rnf) - (n n 2 + n n 2 T n k 2 )] [(n 22 2 + n 22 2 T n k 2 ) + (n 32 2 + n- 2 T n k 2 ) n k 2 ] 



- [(^ 41 2 + n- 2 T n k 2 ) (n u 2 + n u 2 T n k 2 ) + + n~ 2 T n k 2 ) (n^ 2 + n~ 3 2 T n k 2 )} n k 2 . (7b) 

Here, r is defined in Eq. ()10b|) later, Q k 2 = (k/u) 2 , and furthermore, Q^ 2 and £l~ 2 T are 
defined by 



a °' a a °' a 



^ 2 = El^^. ^ = £^7^; (8b) 
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sZ 31 = 


2 

2 / r / ,2 ^ a ' a ' 


D~2 


2 

V"^ n O,a^o,a , 


^32 = 


2 

n O,a^O,a , 

a 


"32,T 


n O,afo,a , 

~~ 2^ r , ,2 ^ T -"' a ' 

a 


fr 2 ■ 

"41 ■ 


a 


fr 2 

"41 ,T 


~ 2^ r , ,2 tT,o,i 

a 


fr 2 - 

"42 " 


= 2^ r 2^' 

a 


"42,T 


= 2^ r ( ;2^ s ' 

a 


O 2 

iz 43 = 


2_^ r , ,2 ^ a - a ' 

a 


— 2 

1 z 43T 


n O,a v O,a 1 


"44 — 


2 / r , >2 ^ a ' a ' 

a 


s 4 44,T 


~~ 2^ r , ,2 V*r,o,o> 

a 



(8c) 



(8d) 

where r 0ia = (1 — f 2 a ) _1 / 2 is the Lorentz factor, and and ip stand for the dephasing factors 
which can be expressed as 

& + 2I/I/! + (i + i/i/l ) T 

01 = 7 , <PT,a = 7 , 

2z/ - vT- a - v{F 
^1 = — 0— ' ^ = — ' 

^ 1 ^ 2 ^ 
02 = -r, <PT = — 7—, 

C2 C20 

■02 = -^ei02, -0T = -^ei0T, 

_ ^ + 2Z/ (l/ + ^aa^ee) , 6^0,5 + [1 + V {l> + V- aa U ee )\ T 

9a,a — 7 1 <PT,a,a ~ 7 j 



_ 2i/ - £i (V + l^aa^ee) , vT- a - {v + ^aa^ee) T 

Va,o - 7 , ^T,a,a - 7 , 

0=-^ = Cl 2 +^ 2 , (9) 

and the abbreviations are 



v = v ee + £ e i, 



eei 



Ci = l + z> e 2 e , C 2 = 1 + 
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6 = Ci - *A 6 = Ci + A 



(10a) 



T = C20T, T = 2J 7o,a, ?s = 3T ,a + T 0ja ; (10b) 

a 

V*. = (IOC) 

and v ee = v ee /u, v ei = v e i/uj, and T ,a = 3T 0ia /r 0ia . In collisionless limits, it follows that in 
Eq. (jlUaj) v, v\ — > 0, and Ci, C2, £1, £2 — ► 1; for infinitesimal thermal spread, in Eq. ()10b|) 
r, T, T s — > 0; and for symmetrically counterstreaming currents of 1*0,1 = — v ^, in Eq. ()10cj) 
"^aa — — 1- It is noted that the second order terms for thermal correction of the form 
of ~ (f , a n- 2 ) 2 have been neglected, as explained in the Appendix. The corresponding 
condition Tq^VL^ 2 <C 1 turns out to be consistent with the results obtained later, as well 
as, with the aforementioned adiabatic condition of \k/u\v t h, a "C 1, where v t h, a ~ ^/T 0>a is 
the electron thermal speed. When assuming the normalized frequencies v ee and v e i to be 
constants, the dispersion Eq. (JJJ) can be expressed as Re(u; 2 , ^ 2 ) +ilm(uj 2 , fi^ 2 ) = 0, which 
contains, in general, ten complex solutions of u>(k), consisting of five pairs of positive and 
negative solutions. In a special case, they may include purely real and/or purely imaginary 
solutions. 

C. Approximate dispersions for specific cases 

In this section, I investigate some specific cases contained in the general result: first for 
To ia = with (1) only electron-electron collisions and (2) only electron- ion collisions, and 
then for (3) T 0>a 7^ and (4) T 0ia = 0, without collisions. 

1. The case including electron- electron collisional effects 

For the case of v ee 7^ 0, v ei — > 0, and Tq^ — > 0, the dephasing factors of Eq. © asymp- 
totically lead to 
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r,2 



, , 1 , 1 + 2 (1 + V aa) V ee 
01, 02 I! 05,a ~> : : . ~ 2 ! 0T, 0T,a, 0T,a,a ~> 0; 

^1, V>2 -> 0; Va,a -» - A ^T-l^ 6 ; V>T, V'T.a, ^T,a,a ~> 0. (11) 

1 + 4z/ ee 

Therefore, in Eq. (|SJ) we read fi^ 2 , fi 22 2 , fij 2 2 , — > 0> an d Eq. © then reduces to 

Re(^ 2 , n- 2 ) « (1 - fi n 2 ) [(1 - fi 2 - 2 ) - (1 + fi 3 - 2 ) fi~ 2 ] - Q^nf, (12a) 

im(^ 2 , fi- 2 ) « - (i - fi n 2 ) ^ fc 2 - n^n^ 2 . (12b) 

Instead of Eq. (JHJ), I introduce the definitions of 



a 



,2 



n /-2 _ \^ ra 0,a^0,a n /-2 _ \^ ^O.a^O.a / 1U \ 



a 

where u' 2 = u 2 / (<p a>a + z'V'aa), which may be rewritten as 



a 

,/2 



^ = l + 2gee ^ (M) 

1 + 2 (1 + Vaa) Vee 

The definitions of fi^ 2 , fi^ 2 , and fij 2 in Eq. (|13a|) are recalled later. It is noted that for 
a trivial case of copropagating currents with u ,i = t>o,2, i-e., v aa = 1, Eq. (JT3j) reduces to 
u/ 2 = c<j 2 , indicating that indeed, dephasing effects vanish. Making use of Eqs. ()13j) and (JHJ), 
the approximate dispersion Eq. (|12j) can be written in the form of 

u 2 (i - ^r 2 ) (i - ^ 2 ) - k 2 [(i - fir 2 ) (i + ^r 2 ) + nfn'i 2 ] « o. (is) 

Equation (|THj) contains six solutions of cu. One should note that the factor of (1 — fij" 2 ), 
regardless of the transformation of Eq. (|14j) . involves a simple eigenmode of the relativistic 
plasma oscillation. More details are discussed in Sec. IIII Al 
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2. The case including electron-ion collisional effects 

For the case of v ee — > 0, i> e i ^ 0, and T 0ja — > 0, the dephasing factors of Eq. (jUJ) asymp- 
totically lead to 



01, 4>a,a -> 02 = -r! 0T, 0T,a, 0T,a,a ~> 0; 
S2 

^a,a ^4>2 = —fr\ ifo, lfo,a, ^T,a,a ~> 0. (16) 
S2 

Therefore, in Eq. (JBJ) 0^ 2 T — > 0, and Eq. (JJJ) then reduces to 



Re(^ 2 , 1\7 2 ) « (1 - fi n 2 ) [(1 - ^2?) - (1 + ^si 2 ) 

— f2 12 (^22 "I" ^32 ^fc ) — (^41 ^43 — ^42 ^44 ) ^fc ' (17a) 

im(^ 2 , n k 2 ) » -q- 2 2 [(1 - fi 21 2 ) - (1 + fi 3 - 2 ) fi~ 2 ] 

- (1 - fi n 2 ) + ^"^fe 2 ) - (^4~4 2 + ^ (17b) 

Here, we note the relations of 



f2^2 ~ VeiQii , f2 2 2 ~ Vev^>2\ 1 ^32 ~ ^ei^31 5 



O" 2 Rj -7/ O" 2 CT 2 Rj -v O" 2 CT 2 Rj O" 2 

ic, 42 t/ ei i4 41 j i4 44 ~ t/ ei i4 43 > i4 42 i4 44 ■ v ±0 / 

For heuristic ways, I properly choose the dephasing boost-frame defined as 

u" 2 = (1 + ^V 2 , (Vw 2 ); fc" 2 = (1 + t9 ei )k 2 } (VA; 2 ). (19) 

Note that in contrast with Eq. ([14)1 . the wave number is also transformed by the operator 
in Eq. ()19j) . Taking account of the transformation, I give the definitions of 
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2 

n //-2 _ n O,a V 0,a Q//-2 _ ^."^O.a ^OfTl 

instead of Eq. (jHJ). By using Eqs. (flT?j) and (j2TJj) . the approximate dispersion Eq. (|T7|) can be 

expressed as 



uj" 2 (1 - fi'r 2 ) (1 - ft 2 '- 2 ) - A;" 2 [(1 - fi'r 2 ) (1 + ft' 3 '- 2 ) + ft 4 '~ 4 ] « 0, (21) 

where f2 4 -4 = (f2 4 ~ 2 ) 2 . Equation (|21j) contains six solutions. It is noted that the form of 
Eq. (|21j) is quite similar to that of Eq. (|15|). though there exist some differences in dephasing 
property between the two. For example, the eigenmode relevant to the plasma oscillation, 
which is contained in the factor (l — fl"~ 2 ), is now undergoing the transformation of Eq. (|T^|) . 
More on these is given later in Sec. IHI Bl and IIIIEI 

3. The case including thermal corrections 

For the case of v ee , v ei — > 0, and T 0)a 7^ 0, the dephasing factors of Eq. Q asymptotically 
lead to 



01, 02, 0a,a ~^ 1} 4>T ~> T = -2T] Tj5 , (j) Tfi ,a ~> ~ [T + T 0( aJ J 

^1,^2, ^a,a, V'T, ^T,a, ^T,a,a ~> 0. (22) 

Therefore, in Eq. (jHJ), fi^ 2 , fi^ 2 T , fij 4 2 , fij 4 2 T — > 0, and Eq. (J7j) then reduces to 

Re(u; 2 , n- 2 ) « [(i + rfi fc - 2 ) - (fi n 2 + ^fV^ 2 )] 
x [(i + rfi, 2 ) (i - n k 2 ) - (n 2 2 + n 21 2 T n k 2 ) - (n n 2 + n^nf) n k 2 ] 

— (f2 41 + f2 41T f2 fc ) (f2 43 + f2 43T f2 fc ) Cl k , (23a) 
Im(cj 2 ,fi fc 2 ) « 0, (23b) 
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and we have the relations of 



In this case, dephasing effects disappear, and the dispersion relation yields purely real and 
purely imaginary solutions. Concerning Eq. ()24|). I give the definitions of 



-L 0,a^ 1 ,a^ 

^ 2 = Ef?^' ^ = rfi 2 ~ 2 ; (25b) 

2 2 
^3 =Z^-F-^' fi 3,T = 2^ 2-^,a,a; (25c) 

n * =Z^-TT^' fi 4,T = 2^ r^7±^,a, (25d) 

instead of Eq. (jHJ). Note that the definitions of fi^ 2 , ^2~ 2 > an d ^4 2 have already appeared 
in Eq. (|13a|) . Using the definitions of Eq. (J25)l . the approximate dispersion Eq. can be 
written as 



[(i + r^ 2 )-(^ 2 + n- 2 T n k 2 )} 
x [(1 + m k 2 ) (1 - n 2 2 - n k 2 ) - (n 3 2 + n^ 2 ) n^ 2 } 



- (Qf + n^ 2 T nf) 2 nf «0. (26) 

In general, Eq. contains ten solutions. At first glance, the form of Eq. ()26|) seems to be 
different from that of Eqs. (|15|) and (|2ip. Indeed, there are some differences in dispersive 
property, for instance, the plasma oscillation mode contained in the factor of [(1 + tQ^ 2 ) — 
{fli 2 + Q^Q^ 2 )} involves the thermal dispersion. In the special case without thermal 
corrections, however, Eq. ()26|) recovers the same form with Eq. (|21|). as shown below. 
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4- The collisionless case without thermal corrections 



For the case of z> ee , v ei — > 0, and To,a — > 0, the dephasing factors of Eq. (0) asymptotically 
lead to 



n, <P2, <Pa,a 



?T, <PT,a, <PT,a,a 



0; 



Ipl, 1p2, 1pa,a, ifo, 1pT,a, 1pT,a,a ~> 0. (27) 

Therefore, in Eq. (JKJ), fi^ 2 , fij 4 2 , fiZ- 2 T — ► and fij 2 « fij 3 2 , and then Eq. (J7J) reduces to 



(i - ^r 2 ) (i - fi 2 ~ 2 ) - fc 2 [(i - (i + fi 3 ~ 2 ) + ^~ 4 ] « o, 



(28) 



where fi^ have been defined in Eq. and fij 4 = (fij ) 2 - It is found that Eq. (J2HJ) 



maintains the form similar to Eqs. 
with that obtained by Califano et al. 




This dispersion Eq. (}2*S|) exactly coincides 



III. EIGENMODE PROPERTIES AND GROWTH RATES OF THE 
RELATIVISTIC CURRENT FILAMENTATION INSTABILITY IN A 
COLLISIONAL PLASMA 

In the following, we find the solutions contained in the approximate dispersion Eqs. (|15|). 
(ED), (|2b|h and ()28|). The complex eigenmodes are explicitly written down, and surveyed for 
wide parameter ranges of counterstreaming relativistic currents. 

A. Eigenmodes including electron-electron collisional effects 



At first, I seek the solutions of Eq. (|15|) in terms of u. Let us consider the symmetrical 



configuration of counterstreaming currents such as no,i = ^0,2 = 0.5 |33|, having t> 0) i = — t>o,2- 
The choice of the parameters may be instructive for making a direct comparison between 
the present results and the previous ones [27(. Equation (JT5J) can be then cast to 



1 

T~ 



J 2 u 2 - 



r 3 
1 



==T I - ( 00' 2 + ^ ] k 



r, 



(29) 
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where v = \v 0ja \ and T = (1 — Vq) 1 / 2 , and the transformation Eq. (fTl| reduces to 



uj 
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(1 + 2iV ee ) uj 2 . (30) 

The first factor of (left-hand side) (lhs) of Eq. (f2U|) yields a simple electrostatic mode, 
corresponding to the relativistic plasma oscillation: uj = ±u r and u>i = 0, where 



cov.(r i 



r 1/2 . (3i) 

Below, we refer to this eigenmode as oscillatory mode (o mode). It turns out that momentum 
exchange between electrons and electrons does not disturb the plasma oscillation. 

The second factor of lhs of Eq. ()29|) contains four solutions. The values of uj' 2 are connected 
with uj 2 through the complex operator in Eq. ()30|). Therefore, in contrast to collisionless 
cases, uj 2 values are of complex, to give 



uj (k ,z/ ee ,r ,/x) = - 



k 2 + ^3 + sgn( y u) v / A(/c 2 , u ee , T ) exp [id(k 2 , z/ ee , T )] 



where sgn(/z = /i ± ) = ±1, and 



A(k 2 ,u ee ,T ) = y/B^(k\i>l,T ) + C 2 (k\i> ee ,T ), 



B(k 2 D 2 T\-h* I 2 ( 2F + 4 ^ £ ~ 1) ,2 , 1 



c(k ,z/ ee ,r ) 



i 



(i + ^i) n 

8z> ee (Tg-1) 

(i + 4*> 2 e ) rg 



A- 2 



■d{k ,z/ ee ,r ) = - tan" 



C(k 2 ,v ee ,T ) 
B(k 2 ,i>l,V ) 



(32) 



(33) 



Note the relation of A + B > and A - 5 > 0. 

Moreover, I rewrite Eq. ()32j) in the form of the polar coordinate: uj = ±\uJo\exp(i6o) , 
where ujq and 9q are purely real numbers. For the signs (±) each, we have two solutions 
with /i = yU + and which are referred to as positive mode (p mode) and negative mode (n 
mode), respectively. The squared magnitude and the polar angle of the complex eigenmodes 
are given by 
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uj (k ,z/ ee ,r ,/i) = -a 



k 2 + -3 + sgn(//)a(& 2 , z/ ee , T ) 

J- n 



+ /5 2 (P,z> ee ,r ), 



(34a) 



Oa(k 2 ,v ee ,r Q ,/i) = -tan 1 



respectively, where 



-sgn(/i)/j(fc 2 ,z/ ee ,r ) 
^ 2 + ft + sgn(/i)a(/c 2 , z> ee , T ) 



(34b) 



a(k ,z/ ee ,r ) 



A(P, I > ee ,ro) + J B(A; 2 ,z> e 2 e ,r ) 



> 0, 



(3(k ,z/ ee ,r j 



A(P,z> ee ,r )-5(A; 2 ,z> e 2 e ,r ) 



> 0. 



(35) 



One should note that Eq. (J34b|) restricts its parameter range to — n/4 < 9 < 0, because the 



argument of right-hand side (rhs) of Eq. (j34b|) is negative definite, i.e., the ratio is — /+ for 
the p mode and +/— for the n mode. In order to reconstruct the polar angles consistent 
with the argument, I appropriately define ujq p = \oj (h = fA+)\ and 9 0p = 6> (/i = (/,+) for 
the p mode, and uj 0n = \uo(/jl = fi-)\ and 9 Qn = #o(/z = //_) + n/2 for the n mode. Note 
the allowable parameter ranges of — 7r/4 < 9q p < and 7r/4 < 9q h < tt/2. For the p 
mode with < sin^op < and V2/2 < cos#op < 1, the complex eigenmodes of cu = 

uj 0p exp(i9 p) and uj = ~u 0p exp(i9 0p ) = uj 0p exp[i(9 0p + n)} reflect decaying and growing waves, 
respectively. Thus, the phase of the decaying and growing waves propagates toward y and 
— y direction, respectively. On the other hand, for the n mode with y/2/2 < sin# n ^ 1 an d 
< cos^on < V2/2, the complex eigenmodes of to = uo n exp(i9 0n ) and uj = —ujQ n ex\>{i9Q n ) = 
a; On exp[z(0 O n — n )] reflect growing and decaying waves, respectively. In contrast to the p 
mode, the phase of the growing and decaying waves propagates toward y and —y direction, 
respectively. 

In Fig.^for the typical current speed of vq = 0.9 (Tq = 2.29), I show the polar coordinate 
plots of the complex eigenmodes u for given wave numbers k, varying electron-electron 
collision parameter v ee . The trajectories of u can be compared to those of the arrowhead of 
vectors. In the collisionless limit of v ee — > 0, we read A w B and C — > in Eq. (J33)) . and 
(3 — > in Eq. (|35p. so that 9 — > in Eq. (|34b[) . Hence, the vectors of the p mode direct the 
real axis, reflecting the purely oscillatory wave mode (9 0p = 0; 9 0p + tc = tt), whereas the 
vectors of the n mode direct the imaginary axis, reflecting the purely growing (9o n = tt/2) 



16 



and purely decaying (0 O n — vr = — 7r/2) mode. The unstable mode just represents the 
electromagnetic Weibel instability in a collisionless plasma, which is discussed in Sec. IIII Dl 
later. 

As shown in Fig. ^a) for k = 3 x 10~ 3 , for the small but finite value of v ef , the vectors of 
the n mode depart clockwise from the imaginary axis. As v ee increases, the real components 
of the vectors increase, while the imaginary components decrease. In this aspect, such a 
growing wave mode is considered to be the dephasing Weibel mode with reduced growth 
rate. In the strong collisional regime, the polar angles of the vectors of the growing and 
decaying mode approach 9$ n = tt/4 and 9o n — ir = Sir/ 4, respectively, and the vectors 
shrink, reducing both their real and imaginary components. The vectors of the p mode do 
not largely depart from the real axis for the small value of k. 

As seen in Fig.^b), for the moderate value of k = 3 x 10 _1 , the trajectories of the n mode 
are similar to those for the smaller k value. As for the vectors of the p mode, I now find the 
clockwise deviation from the real axis. This represents the decaying (6*o p < 0) and growing 
(6o p + 7r < 7i") electromagnetic mode. In the moderate collisional regime, the magnitude of 
the deviation angles tends to be large, though the maximum value is found to be small, 
compared with that for the n mode. In the strong collisional regime, the vectors of the p 
mode are likely to return to the real axis, reducing their magnitude. For comparison, the 
fixed vectors of the o mode are also plotted in the Fig. d(a,b). In these cases, the magnitude 
of the vectors is larger than that of the p and n mode. 

For the p and n mode, the angular frequency of the oscillations can be defined by u r (fi = 
= uopCosOop > and u r (fi = = cuoncos^on > 0, respectively; and the linear growth 
rate can be defined by 7(/i = = —Uo p sm9o p > and 7(/x = = u 0n sm6 0n > 0, 
respectively. They are summarized as follows: 



u r (k ,v ee ,r ,/i) = - * 2u%(k 2 ,v ee ,T ) + 



k 2 



r 3 

1 



sgn(/x)a(/c 2 , v ee , T ) 



(36a) 



7(fc 2 ,z/ ee ,r ,/i) = -\j2u^(k 2 ,i> ee ,T ) - 



k 2 + -3 + sgn(fi)a(k 2 , u ee , T ) 

J- n 



(36b) 



Here, note the relation of u r (fi + ) > 7(/i+) and oj r {jiJ) < 7(/i_). 

In Fig. El for vq = 0.9, I show the growth rate 7 of Eq. (|36b|) as a function of the wave 
number k for given collision parameters v ee = 0.1, 1, and 10. For the n mode, it is found 
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that the collisional effects lower the growth rate for all k, maintaining the dependence of 
7 oc k for small value of k, as well as, the saturation property for large value of k. Such 
asymptotic properties of the grow th rate can be also seen in the purely growing Weibel 
mode in a collisionless plasma 27, 2sj]. In addition, the collisional effects create the p mode 
due to the dephasing mechanism mentioned above. The growth rate has the dependence of 
7 oc A; 2 and 7 oc k~ l in the small- and large-/c region, respectively, taking a peak around the 
moderate value of k. Such a peak tends to be prominent for v ee ~ 0(1), and then decreases 
as v ee further increases, as consistent with Fig. ^b). As a result, the growth rate cannot 
exceed that of the n mode corresponding to the dephasing Weibel mode with reduced growth 
rate. Within the present framework, it seems that both modes do not definitely cut off the 
growth of short wavelength perturbations. 

In Fig. El for v = 0.9 and v ee = 1,1 show the growth rate 7 as a function of the angular 
frequency uo r of Eq. ()36a|) . varying the wave number k as a parameter. The growth rate of 
the n mode turns out to be larger than that of the p mode for all k, which is in contrast with 
the case including only electron-ion collisional effects, as shown in Fig. |H1 later. Moreover, I 

3/9 

found that the angular frequencies of the p and n mode are separated at u r = T Q ~ 0.29, 
and the oscillation frequency of the p mode is always higher than that of the n mode. The 
separation frequency, where duo r /dk —* +0, is lower than the frequency of Eq. ()31|) for the 
mode, i.e. the plasma cutoff frequency. The relation of du) r /dk > for both the p and 
n mode ensures that the direction of the group velocities coincides with that of the phase 
velocities of carrier wave. For large value of k, the oscillation frequency of the p mode 
goes far beyond the relativistic plasma frequency, while the growth rate decreases. These 
properties again appear in the case including electron- ion collisional effects (see also Fig. |HJ). 

In Fig. 13] for vq = 0.9, I show the growth rate 7 as a function of the collision parameter 
v ee for given values of k = 0.01 and 1. It is found that for moderate value of u ee , the growth 
rate of the p mode takes the peak value, which tends to be well pronounced especially for 
k ~ 0(1). For example, for k = 1 the growth rate takes the peak of 7 ~ 5.8 x 10 -2 at 
v ee ~ 0.63. In the weak collisional regime, it has dependence of 7 oc z/ ee , but on the other 
hand that of the n mode is almost constant, to give, e.g., 7 ~ 0.51 for k — 1. In the 
strong collisional regime, the growth rates of both the p and n mode decrease, to exhibit the 
asymptotic behaviors of 7 oc z/" 1 and 7 oc £~ ' 5 , respectively. 
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In Fig. |S] for v ee = 1, I show the growth rate 7 as a function of k for given values of 
v = 0.5, 0.9, and 0.99, corresponding to T = 1.15 (79.0 keV), 2.29 (611 keV), and 7.09 
(3.11 MeV), respectively. For small value of k, both the p and n mode seem to increase their 
growth rates as T increases. For large value of k, however, such properties appear merely in 
weak to mild relativistic regime. Namely, in strong relativistic regime, the growth for large 
k tends to be suppressed, for example, as for the n mode, the saturation level of the growth 
rate decreases as To increases. The peak of the growth rate of the p mode is likely to shift 
to the smaller-Zc region, reducing its value. Although, the energy dependence of the n mode 
appears again in the case including electron-ion collisional effects, the p mode significantly 
changes its property, as shown below. 



B. Eigenmodes including electron-ion collisional effects 

I seek the solutions of Eq. (|2"T|) in terms of u. In the case of the symmetrically counter- 
streaming currents with n ,i = ^0,2 = 0.5 and vqi = —^0,2, Eq. (J2T|) can be cast to 



LU" 2 ' 



112 I 112 1 \ ( 112 , V \ K//2 



, , - , , , , 0, (37) 

where vq = \vo >a \ and r = (1 — Vq) ~ x l 2 . Note that Eq. (|3*7jl has the same form as Eq. (|2T?|) 
for the previous case. However, now all u 2 and k 2 values are being dephased by the transfor- 
mation Eq. ((EH) , to provide the eigenmodes significantly different from those derived from 
Eq. (|29|). The first factor of lhs of Eq. (|37|) contains a modified electrostatic mode of the rel- 
ativistic plasma oscillation. The eigenmode may be written in the form of uj = ±(av — iuii), 
and the growth rate can be then defined by 7 = Ui, that is, 



"r&ei, To) = 3/4 , 1{V«, T ) = " (38) 

vroC 2 iKi) vroC 2 v4o 

It is found that the dissipative effects owing to the electron-ion collisions create the growing 
and decaying oscillatory mode, which seemingly carries its phase toward — y and y direc- 
tion, respectively. As far as ignoring thermal and asymmetrical effects of counterstreaming 
currents is concerned, this mode is independent of wave number, so that the group velocity 
is null. Note that in the collisionless limit of u e i — > 0, Eq. (|3*H|) reduces to Eq. ()31|) which 
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denotes the relativistic plasma oscillation. In this sense, we refer to the eigenmode identified 
by Eq. (J3"8j) . for convenience, as o mode. 

The second factor of lhs of Eq. (J37|) contains four complex solutions. Along the 
manner explained in Sec. 1111 Al I express the solutions in the polar coordinate form of 
uj = ±\u \exp(i9 ) , where u and 9 are purely real numbers. For the signs (±) each, we 
get two solutions. The squared magnitude and the polar angle of the complex solutions are 
given by 



u (k ,u ei ,r Q ,n) = — 



(2k 2 + =3 + sgn(/i) (a + u ei (3) 
1 



+ 



-f + sgn(/i) (u ei a - (3) 
1 



(39a) 



9 (k 2 ,u ei ,T ,fi) = -tan 



1 



sgn(/i) {v ei a - (3) 



(2k 2 + 4 + sgn(/i) (a + u ei (3) 



respectively, where 



(39b) 



a(k 2 ,u ei ,T ) 



A(k 2 ,i) ei ,T ) + B(k 2 ,i) 2 t ,T ) 



> 0, 



and, 



(3(k ,z/ ei ,r ) 



A{k 2 1 p eil Y )-B{k 2 ,ul,Y, 



>0, 



(40) 



A(k 2 , P ei , To) = JB 2 (k 2 , u 2 ei , T ) + C 2 (k 2 , 9 ei , T ), 



B(k 2 , i£, r ) = (1 - v 2 ei ) k* + 2 (2r ; 3 l) k 2 + 1, 

J- n J- n 



C(k 2 ^ eu To) = ^k 2 [k 2 + ^ 



2TI-1 



(41) 

Note the relation of A + B > and A — B > 0. Below, I refer to the solutions with fi = fi + 
and as p mode and n mode, respectively. As is the case with Eq. ()34b|) . Eq. ()39b|) 
holds a parameter range of— 7r/4 < 9 < 0, because the argument of rhs of Eq. (|39b|) is 
negative definite, i.e., the ratio is — /+ for the p mode and +/— for the n mode. Thus, 
recalling the definitions of co> 0p = \uq(/jl = fi + )\ and 9 0p = 9 (p, = fi + ) for the p mode, and 
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uon = \uo(fj, = //_)| and 9o n = = A 4 -) + 7r /2 for the n mode, the complex eigenmodes 
to = ±uj 0p exp(i9o p ) denote the decaying (+) and growing (— ) wave that carry their phases 
toward y and — y direction, respectively; and uj = ±uj 0n exp(i9 0n ) denote the growing (+) 
and decaying (— ) wave that carry their phases toward y and —y direction, respectively. 

In Fig. |U1 for the current speed of Vq = 0.9, I show the polar coordinate plots of the 
complex eigenmodes uj for given wave numbers k, varying electron-ion collision parameter 
v e i. We compare the trajectories of uj to those of the arrowhead of vectors. In the collisionless 
limit of u ei -> 0, we read A « B and C -> in Eq. (JUJ), and ft -> in Eq. (|£0j). so that 
#o — ► in Eq. ()39bj) . Hence, the vectors of the p mode direct the real axis, reflecting the 
purely oscillatory wave mode (8q p = 0; 8 0p + tt = tt), while the vectors of the n mode 
direct the imaginary axis, reflecting the purely growing (8 0n = tt/2) and purely decaying 
{Oon—n = — tt/2) mode. The unstable mode represents the electromagnetic Weibel instability 
in a collisionless plasma. 

As shown in Fig.lH^a) for k = 3 x for small but finite value of v e i the vectors of the p 
mode depart clockwise from the real axis. As v e i increases, the real components of the vectors 
decrease, while the imaginary components increase. The magnitude of the deviation angles 
is larger than that in the case including only electron-electron collisional effects [compare 
Fig. ma)]. In the strong collisional regime, the vectors shrink, reducing both their real 
and imaginary components. This is the decaying and growing electromagnetic mode, which 
possesses the allowed range of polar angle of — tt/4 < 9 0p < and 37r/4 < 9 0p + tt < tt, 
respectively. For the parameter range of v e i being considered, the vectors of the n mode 
do not largely depart from the imaginary axis, that is, the purely growing Weibel mode is 
not so dephased. The most remarkable property can be seen in the o mode. For v ei ^ 0, 
the vectors leave the real axis, and as v e i increases the magnitude of the deviation angles 
increases. As consistent with Eq. ((SB)) , at v ei = 1 the vectors of the growing and decaying 
oscillatory mode have the angles of 37r/4 and — vr/4, respectively; and in the strong collisional 
limit, asymptotically approach tt/2 and —tt/2, respectively. 

In Fig. Efb), for the moderate value of k = 3 x 10 _1 , I now clearly find, for v e i ^ 0, the 
clockwise deviation of vector pairs of p, n, and o modes. This essentially means that all 
modes are in phase lag, because of the frictional nature of collisions. The vectors of the p 
and n mode exhibit the behavior similar to that displayed in Fig. db). In contrast with the 
case for small k, the polar angles of the vectors of the p mode cannot reach to 37r/4 and 
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— 7r/4, and in the strong collisional regime the vectors are likely to return to the real axis, 
reducing their magnitude. On the other hand, the vectors of the n mode largely deviate 
from the imaginary axis until the angles reach to 7r/4 and Sir/A. In the weak collisional 
regime, the vectors increase the real components and decrease the imaginary components, 
whereas in the strong collisional regime, there is decrease both in the real and imaginary 
components. As a result, the electron-ion collisional effects lower the growth rate of the 
n mode, namely, the dephasing Weibel mode. Note that the trajectories of the vectors of 
the o mode are the same as those displayed in Fig. |Ufa). For the moderate value of k, the 
magnitude of the vectors of the n mode is found to be comparable to that of the o mode. 
For the p and n mode, the angular frequency of the oscillations can be defined by u r (fi = 
= uj 0p cos9 0p > and u; r (/i = = u Qn cos9 0n > 0, respectively; and the linear growth 
rate can be defined by 7(/i = = — uj 0p sm9 0p > and -y(fi = = ^onsin^on > 0, 
respectively. They are summarized as follows: 



uJ r (k 2 ,u ei ,T ,fi) = - 



\ 



2u 2 (k 2 ,i) ei ,T ) + 



^ + sgn(» [a(k 2 , u ei , T ) + is e if3(k 2 , v e i, T )] 
k 2 H — - 



(42a) 



"f{k 2 ,v ei ,r ,(i) = - 



\ 



2u 2 (k 2 ,u ei ,T ) - 



k 2 1 



k + sgn(/i) [a(k 2 , P ei , T ) + v e iP(k 2 , v ei , T )] 



(42b) 

Here, note the relation of u r ((/,+) > 7(/i+) and u r (/j,-) < 7(/i-). 

In Fig. [3 for vq = 0.9, I show the growth rate 7 of Eqs. (fHKJ) and (J42bj) as a function of 
the wave number k for given collision parameters v e i = 0.1, 1, and 10. For the n mode, the 
collisional effects lower the growth rate, especially in the large-fc region where the growth 
rate tends to saturate. In contrast to the case including only electron-electron collisions, 
the growth rate for small k, which has the dependence of 7 oc k, is found to be not so 
depressed due to the electron-ion collisional effects (compare Fig. EJ). The p and mode, 
which appear owing to the collisional effects, are prominent for the parameter of v e i ~ 0(1), 
and sufficiently surpass the n mode in growing long-wavelength (small-/c) perturbations. 
This feature is in major contrast to the case including only electron-electron collisional 
effects. For the larger collision parameter, it is noteworthy that even in the large-fc region, 
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the growth rate of the o mode exceeds that of the n mode, i.e., the suppressed Weibel 
mode. This is one of the most important results in the present paper. In the large-fc 
region, the growth rate of the p mode decreases, showing the dependence of 7 oc and 
is far below that of the n and mode. In the small-fc region, the growth rate of the p 
mode is almost independent of k, and is always smaller than that of the mode. That is, 
7 (^,r ,/i = fi+) « [(v^2 - l)/(2( 2 Tl)} 1 / 2 < zV(v/IX 2 3/4 ) for v ei > and T > 1. 

In Fig. |S1 for v = 0.9 and v ei — 1, I show 7 of Eq. (|42b|) as a function of the angular 
frequency uo T of Eq. ([42 aj) . varying k as a parameter. For comparison, in the u r — 7 plane, 
I also plot a fixed point given by Eq. (|H8|) for the mode. It is found that for large value 
of k, the growth rate of the n mode is larger than that of the mode, which is always 
larger than that of the p mode. Note that, for the p and n mode, duo T /dk > 0, and the 
angular frequencies of both modes are clearly separated at the frequency of u r ~ 0.2, where 
du r /dk — > +0. As a result, the oscillation frequency of the p mode is always higher than 
that of the n mode. These features could be also seen in Fig. El for the case including 
electron-electron collisional effects. 

In Fig. El for t> = 0.9, I show the growth rate 7 of Eq. (jSHJ), and (J42bj) as a function of 
the collision parameter v ei for given values of k = 0.01 and 1. In the weak collisional regime, 
the growth rates of the and p mode are both proportional to z/ ei , while the growth rate 
of the n mode is almost constant. For smaller k, the growth rates of the and p mode can 
more readily exceed the growth rate of the n mode, and for the moderate value of u e i, they 
take the peak values. The growth rate of the mode can, even for large k, exceed that of 
the n mode, and takes the peak of 7 ~ 0.41 at v e i ~ 1.4. The growth rate decreases as v e i 
further increases, showing the dependence of 7 oc £~ ' 5 . It is also noted that in the strong 
collisional regime, the growth rates of the p and n mode have the dependencies of 7 oc 
and 7 oc ^°' 5 , respectively. As would be expected, anyhow, all modes are suppressed in the 
strong collisional regime. 

In Fig. for v e i = 1,1 show the growth rate 7 of Eqs. and (J42b)l as a function 
of k for given values of v = 0.5 (r = 1.15), 0.9 (r = 2.29), and 0.99 (r = 7.09). It 
should be remarked that for the lower current speed, the mode becomes most dominant. 
In contrast to the case including only electron-electron collisional effects, the and p mode 
monotonically reduce their growth rates as T increases, and soon the mode is overcome 
by the n mode from large-fc region. The energy dependence of the n mode is similar to 
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that shown in Fig. namely, in weak to mild relativistic regime, the saturation level of the 
growth rate increases as Tq increases, while in strong relativistic regime, it tends to decrease. 
In small-fc region, the growth rate seems to increase as T increases, since its curve, roughly 
proportional to k, shifts to the smaller-fc region. Such an apparent redshift can be also seen 
in the large-/c region for the p mode. 



C. Eigenmodes including thermal corrections 



In Eq. (I26j) . I seek the solutions in terms of uj. Regarding the symmetrically counter- 
streaming currents of n 0) i = n ,2 = 0.5, t> 0) i — — ^0,2, and T ,i — ^0,2, Eq- (J2H1) can be cast 
to 



x 



2 / 2 

UJ UJ 



,,r I uj 2 - AT k ) — — I ,.r 



3T n k 2 



ATnk 2 



UJ 



r 3 
1 



To 



<u 2 



3T k 2 k 



(43) 



where vq = \vo >a \, To = (1 — v o)~ 1 ^ 2 i an d Tq = T 0ta . The first factor of lhs of Eq. (|4^|) contains 
four purely real solutions, but two of them are found to be inconsistent with the assumption 
of TqQ^ 2 <C 1 (not shown). The other two solutions can be expressed as u = ±u r , and for 
Ar T k 2 < 1, we obtain 



uj r (k 2 , Tq, Tq 



1 



-T k : 



where T n 



/r V 2 

T 0>a . Note that Eq. is of the order of uj. 



(44) 

and therefore, the 



assumption TqVL^ 2 <C 1 requires 3T k 2 <C 1, which is consistent with the aforementioned 
relation of ~ r T A; 2 <C 1. Equation (j4^j) just corresponds to the relativisticall y e xtended 
dispersion of the Bohm-Gross wave with nonzero group velocity (e.g., see Ref. |48j for the 
nonrelativistic limit). In contrast to Eqs. (}3T|) and (|3*Hj) . we find the thermal dispersion 
terms characterized by ~ T k 2 in Eq. (|44j) . Physically, this reflects the Debye screening by 
electrons In the limit of T k 2 — > 0, Eq. (j4^|) reduces to Eq. (pITj) for oscillatory mode. 

Moreover, the second factor of lhs of Eq. (|4*3*|) contains six solutions. Numerical calculation 
indicates that these consist of four purely real solutions and a purely imaginary solution 
concomitant with its complex conjugate: uj = rkiuji. For the purely growing mode that is 



24 



of interest here, we can define the linear growth rate as 7 = Ui. In Fig. ^2 (solid curve), I 
show 7 as a function of the wave number k for vq = 0.9 and To = 0.1, as an example. It is 
found that the thermal effects simply lower the growth rate of the purely growing Weibel 
instability, at least, in the range of k ~ 0(1). Note that in contrast to collisional cases, 
the thermal effects do not take part in dephasing the purely oscillatory and purely growing 
mode, but merely give rise to the mode- dispersion which can reduce the growth rate. 



D. Eigenmodes of collisionless case without thermal corrections 



In the collisionless limits without thermal corrections, the dephasing operators in Eqs. (|T9"j) 
and (|30|) asymptotically approach unity, vanishing their imaginary parts. Then, Eqs. 
(J37J), and (jlSj) degenerate into a unique equation, 



to 



r 



2 / 2 

UJ UJ 



r 



0. 



(45) 



where v = \v 0a \ and T = (1 — Vq)~ 1//2 . As expected, Eq. fl43|) has the same form with 
Eqs. (|2*H|) and (jHZj). The first factor of lhs of Eq. ()45)1 yields the plasma oscillation mode 
given by Eq. (p?Tj) . 

The second factor of lhs of Eq. (|43j) contains four solutions. In contrast to collisional 
cases, the values of uj 2 are obtained as purely real numbers. Concerning the polar form 
uj = ±|c<Jo|exp(i#o), Eqs- p4f and (}39|) reduce to 



uj (k ,r ,/i) = - 



fc 2 + 4 + sgn( /U ) v / A(P,r c 



r 3 
1 



(46a) 



9 = 0, (46b) 

where, 

A(k\r ) = k*+ H2T °- l) e + ± (47) 

1 1 

Recalling the definitions of uo p = |^o(a* = / i +)l an d 0o P = #0 for the p mode, and u Qn = 
\oJo(n = /i_)| and 9o n = 9 + n/2 for the n mode, Eq. ()46bj) leads to 9 0p = and 9 0n = tt/2. 
It is, therefore, found that for the p mode, the real solutions of u = ±u 0p exp(i9 0p ) = ±ojq p 
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describe the purely oscillatory wave mode, while for the n mode, the imaginary solutions 
of u — ±uo n exp(i9 0n ) = ±iuo n describe the purely growing (+) and purely decaying (— ) 
mode. The phase properties are shown in Figs. Q and El 

Now, the angular frequency of the oscillation and the linear growth rate can be defined 
by uj r = Lu p > and 7 = uj 0n > 0, respectively. Taking the inequality y/~A > k 2 + (1 /Tq) 
into consideration, the angular frequency and the growth rate can be expressed as 



uj r (k 2 ,r ) = — 



v / I(F 2 ,r )+(F + f 3 



l(k 2 ,T c 



1 

71 



k 



r 3 
1 



1/2 



1/2 



(48a) 
(48b) 



VA(k 2 ,T )- 

respectively. Note the relation of uj r > 7. In Fig. ^2 (dotted curve), I show the linear growth 
rate 7 of Eq. (J48bj) for vq = 0.9, as a function of the wave number k. This growth rate is just 
of the relativistically extended electromagnetic Weibel instability in a collisionless plasma 
It is found that for T^k 2 -C 1 and T^k 2 ^> 1, Eq. (j48b|) simply exhibits the asymptotic 



3- 



property of 7 ~ a/T^ — lk and 7 ~ a/ (Tq — 1)/Fq, respectively. It might be instructive to 
compare Fig. ^2 with Figs. El and for the case including collisional effects. It is confirmed 
that this mode is in a special case of the modes presented in Sec. IIII AI - ITTlCI 



E. Eigenmodes including asymmetric effects of counterstreaming currents 

For another comparison, we are concerned with an asymmetrical configuration of coun- 
terstreaming currents. First, in Eq. (J28| for the collisionless case, I change the parameters 
to an asymmetrical, though still current-neutral initial beam configuration with no,i = 0.1 
(^0,1 = 0.9) and n ,2 = 0.9 (^0,2 = —0.1). Note that the corresponding beam-to-plasma 
density ration, rii,/n p ~ 10, is of relevance to electron transport in the context of ignitor 
physics. Here the electron beam density could be estimated as ~ n c w 10 21 cm -3 and 
the plasma density as n p w 10n c |37| . It is around the low-density foot of steeply rising 
density profile of laser-ablative plasma j^], where the filamentation dynamics may be most 
prominent. As shown in Fig. El (crosses), it is found that for such parameters the growth 
rate reduces by a factor 10 in the small-A; region and by factor vT(J in the saturation region. 
The strong reduction is, more or less, favorable for energetic electrons to propagate through 
the ablative corona surrounding a highly compressed ignitor plasma. 
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In the denser region, electron-ion collisions might play an important role in attenuating 
slow return currents, owing to the relation of [yuj '^ v b){yeii ^pe)" 1 — {nb/n e ) 5 ^ 2 T b 3 ^ 2 <C 1, 
where and lv p i, denote the beam electron-ion collision frequency and beam electron plasma 
frequency, respectively. At this juncture, in order to take account of the collisional effects 
more plausibly, one may replace the dispersion Eq. (J21|) with 



uj 2 



(i - irr 2 ) (i - ir 2 '- 2 ) - k 2 [(i - ft?- 2 ) (i + ft' 3 '- 2 ) + ^r 4 ] « o, (49) 

where, instead of Eqs. (fT§|) and (|2Uj). use u;" 2 = (1 + iu eita )uj 2 and 



respectively. For simplicity, we set the current parameters to the same as the previous ones: 
rio,i = 0.1 (i>o,i = 0.9) and no,2 = 0.9 (t>o,2 = —0.1), but now providing the slow return current 
with r ,2 — 1 to be resistive such as v ei ^ ^ (3> z^o)- Numerical calculations indicate that 
Eq. (|49|) contains six complex solutions, consisting of three pairs of the positive and negative 
solutions. For the three unstable modes each, one can define the growth rates as 7 = \uji\. 

In Fig. I show 7 as a function of k, for the collision parameters of z/ ei l = and 
v e i,2 = 1, as an example. It is found that as a whole the wave number dependence itself 
is similar to that for the symmetrically counterstreaming cases (compare Figs. 171 and fTUj) . 
The purely growing mode {v e i,2 — 0: crosses) is disturbed due to the collisional effects on 
the slow return current to yield the n mode, namely, the dephasing Weibel mode. However, 
the dephasing effects result in only a slight increase of the growth rate. As shown in Fig. El 
(solid curve) for z/ e j 2 = 1 the increasing rate is about 14% at most around k ~ 0.5. I mention 
that for the larger collision parameter, the growth rate of the n mode is depressed below that 
for the collisionless {v e i,2 = 0) case. In addition, the dephasing effects create the unstable 
modes analogous to the and p mode that were introduced in Sec. 1111 Bl The growth rate 
of the corresponding mode now quite weakly depends on the wave number. Both the 
modes are found to surpass the n mode in growing long-wavelength perturbations, as shown 
in the figure. In particular, for the parameter region of ~ 0(1), the growth rate of the 
mode exceeds, even in the short-wavelength region, that of the n mode, that is, along the 

27 



mechanism elucidated in Sec. IIII Bl the electron-ion collision affects the slow return current, 
even if the forward beam current is in the collisionless regime. This consequence is consistent 
with the previous results obtained by carrying out the fully relativistic and electromagnetic 
particle simulation for asymmetrically counterstreaming electron currents in plasma 



IV. CONCLUDING REMARKS 



In conclusion, I have systematically investigated the details of collisional and thermal 
effects on the relativistic current filamentation instability, generalizing dispersion relation 
of the beam-Maxwell system. For specific cases, the approximate dispersions have been 
derived and applied to the instability analysis of typical counterstreaming relativistic 
electron currents, relevant to ignitor physics. For the symmetrically counterstreaming cases, 
the particular results are summarized as follows: 



(i) Effects of electron-electron collision suppress the relativistic Weibel instability for all 
wavelengths. The effects newly create a growing wave mode, but its growth rate is always 
lower than that of the suppressed Weibel instability. 

(ii) Effects of electron-ion collision suppress the relativistic Weibel instability, especially for 
short-wavelength perturbations. The effects create growing oscillatory and wave mode. For 
long-wavelength perturbations, the growth rates of both modes tend to exceed the growth 
rate of the suppressed Weibel instability. For the stronger collisional coupling, and lower 
current speed as well, the growth rate of the oscillatory mode can, for all wavelengths, exceed 
that of the suppressed Weibel instability. 

(iii) Effects of thermal spread simply suppress the relativistic Weibel instability, at least, 
in a moderate wavelength range. 



For the asymmetrically counterstreaming case, the relativistic Weibel instability for the 
symmetrical case is strongly suppressed, though the electron-ion collision still affects the 
slow return current, creating the unstable modes, as mentioned above in result (ii). 

The important point is that, in general, the growing oscillatory mode is created by de- 
phasing a purely oscillatory mode, and the growing wave mode is created by dephasing either 
a purely oscillatory wave mode or a purely growing mode. While the collisional effects invoke 
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phase lag, reflecting the inverse transformation of Eqs. ()14j) and (|19j) for electron-electron and 
electron-ion collision, respectively, the thermal effects do not dephase the purely oscillatory, 
oscillatory wave, and growing mode, but involve mode dispersion. Hopefully, intense laser- 
plasma interaction experiment will be able to reproduce these fundamental consequences, 
though they were derived by leaving out longitudinal modes, complexities of mode-coupling, 
and so forth. 
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APPENDIX A: DERIVATION OF THE GENERALIZED DISPERSION RELA- 
TION INCLUDING COLLISIONAL EFFECTS AND THERMAL CORRECTIONS 

In this appendix, I briefly explain the derivation of Eq. (JJJ). By linearizing the 
continuity Eq. (J2J), the density perturbation of electron component a is described as 
n i,a = n o,a^k \Pi,a,y/To, a ), where Q^ 1 = k/u, and pi >a ,i stand for the first order quan- 
tities of i-directional momentum of the component a. Making use of the relations of 
Pi, a ,x = Ei,x/ ~ v ei)i which are obtained by linearizing x component of the vector Eq. (J3J), 
the first order momenta Pi a ,y can De expressed as functions of the first order electric fields 
Ei^. Substituting these expressions into the linearized continuity equation mentioned above, 
I obtain the first order equations for density perturbation in the form of 



ni t g 
no,a 



Q- 1 v oA l 


(l - f - a Qf + 


+ WeeVaa 


Ei, x 




El,y 


lujT ^ (i - f 0A nf + iu) 


(1 - f ora nf + iu 


+ "1 



(Al) 



for uj ^ 0. Here, the abbreviations of Eq. (fTUj) have been used. Substituting the expressions of 
Pi, a ,i, Eq. (|Aljl . and B\, z = —^Ei^ derived from Eq. (jlj), into the linearized Eq. (J5J), yields 
the first order equations of E^. For the manipulations, the second order terms for thermal 
correction of the form of ~ (T a Q^ 2 ) 2 in the products of, e.g., (1 — T 01 O^ 2 )(l — T 2 ^fc 2 ) are 
neglected, such that (1 - T 0i i^ 2 )(l - T , 2 O fc 2 ) « 1 - (f ,i + f 0>2 )Qf = 1 - Ttt k 2 . Finally, 
we arrive at the equation of the form of Vj l (k,uj)Ei ti rs (VEx^), where the determinant of 
the dielectric tensor is 
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\v(k,u)\ = {(i + rnf) (i - nf) - (q 21 2 + ^ 22 2 ) 

— [(0 21T + i^22,r) + (^31 + ^32 )] ^fc — (^31,T + *^32,t) ^fc } 

^ (1 + rl\T 2 ) - (^ri 2 + - fogr + ggr) "fc 2 

[(^43 + ^44 ) + (^43,T + Z ^44^t) ^fc ] 

- [(n^ + ^ 4 - 2 2 ) + (q 4 - 2 t + zn 4 - 2 T ) nf] n- 1 . (A2) 

Here, the definitions of Eq. (jSJ) have been used. For 7^ 0, the dispersion relation can be 
defined by \U(k,u)\ = 0, to give Eq. ((7|). 

As it is well known, if the imaginary part of the complex eigenvalues u is much smaller 



than the real part, one can calculate them by the Taylor expansion method [5l|. However, 
this is not the case being considered, as seen in Figs. [U El El an d El That is why, I have 
attempted to directly solve the complex Eq. (jJJ) to extract the eigenvalues. 
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FIG. 1: Polar coordinate plots of complex eigenmodes contained in dispersion Eq. (|29J). in the 
parameter range of 1CT 2 < v ee < 10 2 for (a) k = 3.0 x 10~ 3 and (b) k = 3.0 x 10 _1 : u = 
±LUo p ex.p(i9o p ) for p mode (solid curves with filled circles) and uj = ±tvo n exp(i8o n ) for n mode 
(solid curves with filled squares). The p and n mode refer to the definition of, in Eq. (|34ajl . 
ujQp = \u>o(fi = fi + )\ and ujQ n = |wo(m = respectively, and in Eq. (|34b|) . 8q p = 0o(m = M+) an£ i 

—1/2 

#0n = ^o(^ = M-) + 7r /2, respectively. For comparison, I also plot the fixed points of co = ±r o 
for o mode (open circles). The horizontal and vertical line crossing at the center correspond to 
the real and imaginary axis, respectively. The plots of u) can be compared to, e.g., as indicated 
by arrows labeled as u, the trajectories of arrowhead of the vectors, whose magnitude is scaled 
logarithmically by the left axis. Here, I have chosen the parameter of t>o = 0.9, setting u e i = and 
To, a = 0. Note that for (a) k = 3.0 x 10 -3 , the vectors of the p mode do not largely deviate from 
the real axis. For an explanation, see the text. 
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FIG. 2: The linear growth rate 7 of Eq. (|36b|) including the effects of electron-electron collision, 
as a function of the wave number k, for v ee = 0.1 (solid curves), 1 (dot-dashed curves), and 10 
(dotted curves). The p and n mode refer to the unstable modes for /j, = (/,+ and in Eq. (|36b|) . 
respectively. Here, I have chosen the parameter of vq = 0.9, setting z> e j = and To >a = 0. 



35 



1 0° 



0) 

3*1 0" 1 

<D 



o 

o 



1 0" 



1 0" 



10- 4 

1 



n-mode 





/c=1 0" 




o-mode 
(plasma osc.) 



1 o- 2 1 o- 1 10° 
Angular frequency co./co pe 



10 1 



FIG. 3: The linear growth rate 7 of Eq. (|36b|) vs the angular frequency uj r of Eq. (|36aj) for the 
p mode (solid curve) and the n mode (dotted curve), varying the wave number k as a parameter. 
Here, I have chosen the parameters of vo = 0.9 and u ee = 1, setting D e i = and To,a = 0. For 
comparison, the angular frequency of the mode given by Eq. (|31|) is also indicated by a bold 
arrow. 
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Electron-electron collision parameter v 



FIG. 4: The linear growth rate 7 for the p and n mode given by Eq. (|36b|) as a function of u ee 
for k = 0.01 (solid curves) and 1 (dotted curves). Here, I have chosen the parameter of vq = 0.9, 
setting v e i = and Tb ja = 0. 
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FIG. 5: The linear growth rate 7 for the p and n mode given by Eq. (|36b|) as a function of the 
wave number k for vo = 0.5 (solid curve), 0.9 (dot-dashed curves), and 0.99 (dotted curve). Here, 
I have chosen the parameter of v ee = 1, setting i> e j = and Tq a = 0. Note that the dot-dashed 
curves for vq = 0.9 are the same as those shown in Fig. |51 
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FIG. 6: Polar coordinate plots of complex eigenmodes contained in dispersion Eq. Q37|). in the 
parameter range of 10~ 2 < v e i < 10 2 for (a) k = 3.0 x 10~ 3 and (b) k = 3.0 x 10 _1 : u = 
±LUo p ex.p(i9o p ) for p mode (solid curves with filled circles) and uj = ±tvo n exp(i8o n ) for n mode 
(solid curves with filled squares), and Eq. (|38|) for o mode independent of k (hair solid curves 
with open circles). The p and n mode refer to Eq. ()39() for fj, = /i + and respectively, and 
the definitions of ujq p , ujQ n , 8q p , and 6o n are the same as those in Fig. ^ caption. The horizontal 
and vertical line correspond to the real and imaginary axis, respectively. The plots of u can be 
compared to the trajectories of arrowhead of the vectors, whose magnitude is scaled by the left 
axis. Here, I have chosen the parameter of vo = 0.9, setting D ee = and Tq^ = 0. Note that for 
(a) k = 3.0 x 10~ 3 , the vectors of the n mode do not largely deviate from the imaginary axis for 
the parameter range being considered (see text). 
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FIG. 7: The linear growth rate 7 of Eq. ()42bl) including the effects of electron-ion collision, as 
a function of the wave number k, for v e { = 0.1 (solid curves), 1 (dot-dashed curves), and 10 
(dotted curves). The p and n mode refer to the unstable modes for [i = fi + and /i_ in Eq. ()42b|) . 
respectively. Hair lines for each z> e j show the linear growth rate 7 of Eq. (J3SJ) for the o mode. Here, 
I have chosen the parameters of vq = 0.9, setting D ee = and Tq a = 0. 
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FIG. 8: The linear growth rate 7 of Eq. (|42b|) vs the angular frequency uj r of Eq. (|42a|) for the 
p mode (solid curve) and the n mode (dotted curve), varying the wave number k as a parameter. 
An open circle indicates the point given by Eq. (|38|) for the mode. Here, I have chosen the 
parameters of vq = 0.9 and v e i = 1, setting v ee = and To ia = 0. 
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FIG. 9: The linear growth rate 7 for the p and n mode given by Eq. (|42bl) as a function of v e i for 
k = 0.01 (solid curves) and 1 (dotted curves). Hair solid curve shows the linear growth rate 7 of 
Eq. (|38f) for the mode, which is independent of k. Here, I have chosen the parameter of vq = 0.9, 
setting v ee = and Tq^ = 0. 
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FIG. 10: The linear growth rate 7 for the p and n mode given by Eq. ()42b|) as a function of the 
wave number k for vq = 0.5 (solid curve), 0.9 (dot-dashed curves), and 0.99 (dotted curve). Hair 
lines for each vq show 7 of Eq. (|38|) for the mode. Here, I have chosen the parameter of v e i = 1, 
setting D ee = and Tq a = 0. Note that the dot-dashed curves/line for vq = 0.9 are the same as 
those shown in Fig. [7J 
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FIG. 11: The linear growth rate 7 as a function of the wave number k for To = 0.1 (solid curve) 
and To = [Eq. (|48b|) : dotted curve] |27|. Here, I have chosen the parameters of i>o = 0.9, D ee = 0, 
and v e i = 0. 
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FIG. 12: The linear growth rate 7 as a function of the wave number k for asymmetrically coun- 
terstreaming currents of vo,i = 0.9 and vq^ = —0.1 with z/ e i,2 = (crosses) and Ve.i,i = 1- For 
the latter case (£^,2 = 1), solid, dashed, and dot-dashed curves show 7 of the n, o, and p modes, 
respectively. Here, I have set v e i t i = for both the cases. For comparison, I also plot 7 for a 
symmetrical case of vq = 0.9 without collisions (hair dotted curve), which is the same as the dotted 
curve shown in Fig. 1111 
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